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Abstract. We study the Dirichlet problem for first order hyperbolic quasi-linear functional PDEs on a 
simply connected bounded domain of R 2 , where the domain has an interior outflow set and a mere inflow 
boundary. While the question of existence of a solution has already been answered in its predecessor, the 
present paper discusses the uniqueness and continuous dependence on the coefficients of the PDE. Under 
the assumption that the functional dependence is causal, we are able to derive a contraction principle 
C _~) which is the key to proof uniqueness and continuous dependence. Such a causal functional dependence 

appears, e.g., in transport based image inpainting. 
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1. Introduction. The subject of this paper is hyperbolic quasi-linear functional 
PDEs of first order in two space dimensions. We consider the Dirichlet problem 
on simply connected bounded domains Cl and bounded functions uq of bounded 
variation as boundary data. The quasi-linear PDE is stated in the space BV - the 
functions of bounded variation - and formulated in non-conservative form as 

g (c[u](x),Du) = f[u]{x)-C 1 in n\E , u\ dn = u , (1.1) 

where ( . , . ) denotes the scalar product in IR 2 , Du is the derivative measure of u 6 BV 
' and £ 2 the Lebesgue measure on R 2 . The dependence on u of the transport field c 

f--*) and the right-hand side / is thereby of a general functional type. Because the PDE 

— , (1.1) is hyperbolic, we have to rule out the case of characteristic points. For this 

y—i purpose, our requirement is, as in [ ], the existence of a time function T : Cl — > ]R 

which is a global Lyapunov function for the transport field c, i.e., c is assumed to 
00 satisfy the causality condition 

o 

(c[u]{x),VT(x)) > B-\VT(x)\ , B>0. (1.2) 

In addition, T is zero on the boundary and increases towards its maximal level E, 
! '"J the stop set. Consequently, all boundary points are inflow points and E is an interior 

rS outflow set. 

Our interest is the uniqueness of the solution u and its continuous dependence 
on the coefficients c and / of PDE (1.1). In [ , section 6] we have already tackled 
such quasi-linear equations in following way: fix the functional argument of the co- 
efficients by some function v to obtain a linear problem, and the linear theory yields 
a unique solution U[v] depending on v. Now, every fixed point u = U[u] solves 
the quasi-linear equation. We have concluded the existence of a solution as a conse- 
quence of the Schauder fixed point theorem, but we have also given an example of 
non-uniqueness. The latter example demonstrates that we need further assumptions 
on the functional argument in order to obtain uniqueness. 
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Our additional requirement is functional causality. With the time function T we 
have a notion of time on the domain O. We call the coefficients c and f functionally 
causal (w.r.t. time T) if the values c[v] (x) and f[v] (x) depend only on those values 
v{z) which the function v takes on the subset {zefl: T(z) < T(x)} - the past of x. 

Functional causality together with Lipschitz continuity of the coefficients (w.r.t. 
v) carry over to the solution operator U and allow us to derive a contraction principle 
for U. The contraction principle then is the key to show uniqueness and continuous 
dependence of the fixed point solution u = U[u]. 

This strategy has been inspired by Kamont's book [ ]. In [ ] the author considers 
the Cauchy problem for hyperbolic functional PDEs with the functional dependence 
being causal w.r.t. the physical time and proofs the local existence of a unique classi- 
cal solution by employing the Banach fixed point theorem. 

Before we get to the point of the theory we should take a moment and think of 
a possible application. In our case that means transport based image inpainting. The 
term refers to the retouching of undesired or damaged portions of an image. For this 
purpose, we have developed, in [ ], a quasi-linear model for which we can give the 
following rationale: imagine a restorer doing brush strokes in the damaged area O. 
Assuming on the one hand that he only uses color given by the data Uq on dD, and 
on the other that brush strokes go along trajectories x{t) - of a vector field c - which 
constantly carry a single color, we end up with the dynamical system 

x' = c(x) , x(0) = Xn G 30 , 

(1.3) 

u = , w(0) = «o(>o) / 

which describes exactly the characteristics of a linear advection problem with right- 
hand side equal to zero. Ideally, in order to obtain an aesthetic inpainting, the vector 
field c would need to reflect the full expertise of our restorer, which clearly is im- 
possible. But at least the vector field c should be adapted to and hence depend on 
the image u. This consideration aims for a quasi-linear model of type (1.1). For our 
image inpainting model, the dependence of c[u] (x) on u is based on an estimation 
of the tangent vector, which is tangent to the level line of u going through x. This is 
because the brush strokes are supposed to continue level lines of uq which have been 
interrupted by O. Moreover, the estimation of the tangent only depends on already 
known information which causes the transport field c to be functionally causal. 

In [ ] we discussed the modelling and the construction of a fast algorithm, but 
we left open the question for the well-posedness of our inpainting model. This ques- 
tion is answered positively in [ ] as an application of the theory presented in this 
paper. 

Outline of the Paper. In section 2 we review the linear theory and summarize the 
results of [ ]. All the requirements of the linear problem are restated here because 
they will be reused in the later sections. The reader already familiar with [ ] can skip 
this section. In section 3 we take up the quasi-linear problem. We collect the require- 
ments on the coefficients of (1.1) and summarize the previous results of [ ] for the 
quasi-linear case. Section 4 is about the uniqueness. After having established the 
contraction principle mentioned above, we conclude as consequence the uniqueness 
in theorem 4.5. The subject of section 5 is the continuous dependence on the coeffi- 
cients of (1.1). In the proof of theorem 5.2 the contraction principle is again the key 
to the continuous dependence. 
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2. Review of the Linear Case. We summarize the theory of the linear problem 
from [ ], because the quasi-linear theory of the later sections is based on it. We begin 
with the requirements on the data of the linear problem: 

(c(x),Du) =f(x) -C 1 in Q\E, 

"|an = "0/ (2-1) 
(c(x),VT(x)) > p- \VT(x)\ . 

The first set of requirements is concerned with the domain O, the stop set E and the 
time function T, while the second set of requirements is concerned with the transport 
field c, the right-hand side / and the boundary data Uq. 

Requirement 2.1. (Domains) Domains CI c IR 2 are assumed to be open, bounded 
and simply connected and to have C 1 boundary. 

Because of requirement 2.1 the boundary 3fi of a domain is a simple closed C 1 
curve. By 7 : ]R — » dCl we denote a generic regular and periodic parametrization of 
dd. That means 7 G C 1 (1R, 3(1) is surjective and 7'(s) 7^ Vs £ R. Furthermore, by 
I = [s,ii[cR we denote an interval such that 7I j is a generator of 7. 

For our problem time functions are global Lyapunov functions whose range cor- 
responds to a finite time interval. As time is usually a positive quantity which in- 
creases, we define the stop set to be the maximal level of T while in literature the 
stop set is often the minimal level, see e.g. [2]. 

Requirement 2.2. (Time functions) Time functions are of type T : O — > R. The 
upper level-sets of T are denoted by 

Xt>a :={xen: T(x) > A} . 

We assume that time functions T satisfy the following conditions: 

1. T G c(fi). 

2. The boundary ofCl is the start level: T|g n =0. 

3. T incorporates a stop set E with stop time 1; 

a) T(x) x e n\E. 

b) T|e = 1, i.e., E is the maximal level ofT. 

4. T increases strictly from dd towards E. That means that any upper level-set Xt>\ zs 
simply connected and 

Xt>a=X^X VAe[0,l[. 

5. Any proper upper level-set is a future domain: for every A e [0, 1[ the set Xt>a satisfies 
requirement 2.1. Furthermore, the field of interior unit normals to the X-levels 

XT=A=dXT>\, Ae[0,l[ 

ofT is denoted by N : D\E — >■ S 1 . N is required to be continuously differ entiable and 
extendable onto dCi, i.e., N 6 C 1 (fi\E). 
6.* T e C 2 (H), with VT(x) = x e E. 

Remark on 6.*: this assumption is in order to ease things in the passages that 
follow. Because of 6.*, we obtain a simple description of the field N on Q\E by 
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N(x) = VT(x)/|VT(x)| and N is continuously differentiable and extendable onto 

an. 

In part 3 of requirement 2.2 we have assumed that T features a stop set E. Here 
we state the geometric properties of allowed stop sets. 

Requirement 2.3. (Stop sets) We assume that the stop sets E satisfy the following 
conditions: 

1. E is a closed subset of fl 

2. £ is either an isolated point, or a connected set with tree-like structure (no loops). 

3. If E is not an isolated point, it consists of finitely many rectifiable C 1 arcs Ej-.' 

E = U • 

fc=l 

X/ie collection {^k}k=i,...,n * s assumed to be minimal in the number n of arcs, so E is 
decomposed by breaking it up at corners and branching points. 

Furthermore, we require for each arc Ej- that its relative interior Ej. has a given orientation 
by a continuous unit normal n^ : Ej. — > S . 

In the case in which E is not just an isolated point, we also need good behavior 
of the maps T and N close to and on the stop set E. For this purpose, we use the 
following concept of one-sided limits towards z 6 E^: 

If x E B,-(z)\E and if r > is sufficiently small, the projection p of the point x is 
unique. In view of this feature we say a point x 6 B r (z)\E is on the right-hand side 
or plus side (respectively, on the left-hand or minus side) of E^ if 

]x~=Y\ = = ■ (Z2) 

Therewith, a sequence (x n )nehl/ x n £ H\E, tends to z 6 Ej. coming from the plus 
side (respectively, the minus side), in symbols 

x„ z+ (x„ -^z_) , (2.3) 

if the sequence converges towards z and almost all elements x n are on the plus side 
(respectively, minus side). 

With the concept of one-sided limits, the good behavior of the maps T and N is 
summarized in the following requirement. 

Requirement 2.4. (Good behavior at Ej 
2. Requirements on T: 

Let y G E and h 6 S 1 . Lei p = p(j//^) ^ e ^ e ^ es ^ possible order for the asymptotic 
formula 

T(y + rh) = l-0(rP) , r -> 0+ . 

We require that there is a bound q such that sup sup p(y, h) < q. 

yes \h\=l 

2. Requirements on N: 



DIRICHLET PROBLEM FOR QUASI-LINEAR HYPERBOLIC PDES 5 

a) N has one-sided extensions onto the relatively open components ~L k and those exten- 
sions are given by ±n fc : 

N+(y) := Km N(x) , N+(y) = -n k (y) , 

N-(y) : Hm N(x) , AT (y) = n k (y) 

for every y e t k . 

b) The derivative DN has one-sided extensions onto the relatively open components i k , 
i.e., 

(DN) + (y) := x lim DN(x) , (DN)"(y) := lim DN(x) 

exist for every y e Ej.. 

cj |DN| 6 L x (n), i.e., poles of\DN\ at corner-, branching- and terminal nodes o/E are 
integrable. This feature is assumed to hold in the case in which E is an isolated point 
as well. 

So far we have the requirements on domains, stop sets, and time functions. Now 
we turn to the assumptions on transport fields. Here, for a given time function T 
with stop set E, the causality of the transport field w.r.t. T and its good behavior 
close to E are the main concern. 

Requirement 2.5. (Transport fields) Transport fields of type c : Q\E — > R 2 are 
required to satisfy: 

1. c E C 1 (D\E) 2 and c features the following properties: 

a) c and Dc are continuously extendable onto dD,. 

b) IfYL is not just an isolated point, then c and Dc have one-sided limits on the relatively 
open C 1 arcs Ej. o/E: 

c + (y) = lim c(x) and c~ (y) = lim c(x) , 

(Dc) + (y) = lim Dc(x) and (Dc)"(y) = lim Dc(x) , 

for every y e E fc . 

2. (im'£ speed and causality condition: 

a) |c| = 1. 

b) There is a lower bound f> > swc/z that 

jS < (c(x),N(x)> < 1 Vxen\E. 
cj Conditions a) and b) hold for the one-sided limits as well, i.e., 

|c + (y)| = |c-(y)|=l 

and 

/5<(c+(y),N+(y)) < 1 , /5 < ( C -(y),N-(y)> < 1 , 
whenever y belongs to some Ej.. 
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3. Let Zk,k G {1, . . ., m} denote the terminal-, branching- and kink nodes of E. For ecery 
e > s«c/z r7zaf eac/z disk B e (zfc) is compactly contained in O, we define the set 

m 

V E := E U |J B^) . 

fc=l 

a) For every admissible e > 0, i/zere z's a bound M e such that 

\Dc(x)\ <m e , Vie n\y £ 

b) \Dc\ EL 1 (Cl),i.e.,polesof\Dc\atz^,k& {1, . . . ,n} are integrable. 
Finally we need a right-hand side and boundary data. 

Requirement 2.6. (Right-hand sides and boundary data) We assume right-hand sides 
f : Q — > M and Dirichlet boundary data uq : dd — >■ IR fo be / G C 1 (O) and «o G BV (3D). 

Remark: with i/q G BV(3n) we mean that for every regular periodic parametri- 
zation 7 of dCl, the pull-back 7* Mo = u 7 is a periodic BV function. Moreover, 
because 3(1 is one-dimensional, J*uq is BV function of one variable. Fherefor the 
boundary data is essentially bounded. 

In the following sections we assume, if not explicitly stated otherwise, that do- 
mains O, stop sets E, time functions F, and boundary data uq satisfy the require- 
ments above. For the quasi-linear case we will add up the list of requirements on 
transport fields c and right-hand sides / later on. 

Results of [ ]. Problem (2.1) has a unique solution u in BV(Ci) which depends 
continuously on all the data of the problem. Employing the method of characteristics 
the solution is given by 

u(x) = u (tj(T (x),x))+ J f oTj(r,x)dr . (2.4) 



Here, n denotes the backward characteristics which are the solution of 

rj'{.,x) = — Co orj{.,x) , J/(0,x) = x G Q\E . (2.5) 

Fhe forward characteristics - denoted by £ - solve 

£'(.,s) =c o£(.,s), £(0,s) = 7 (s) G 3fl . (2.6) 

According to [ , section 2], Tq = 1 — (1 — F(x)) 1 ^ denotes the transformed version 
of F whose gradient VFd blows up at E. Fo is equivalent to F for it has the same 
level lines as F (level lines get only new names). With Fo we also transform the 
transport field and the right-hand side by c$ = c/ (c, VFo) and fo = f I (c, VFo). 
The forward characteristics w.r.t. cq have then the useful property Fo(£(£, s)) = t, 
i.e., Fq denotes the time of the characteristics. Moreover, Fq features the properties 
Fo|e = l/Folan — 0, and |VFo(x)| >m$> which imply the bound |cq| < l/rwo/3, a 
time range of t G [0, 1] for f ( .,s) , and a bound of l/m^B on the arc length of .,s). 
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By v{t, s) = u o £(f, s) the solution m is represented in characteristic variables 

(f,s) as 

t 

v(t,s) = y*u (s) + J /o°£(T,s)dT. (2.7) 


Theorem 4.1 of [ ] shows that u £ BV(fi) and not only m 6 BV(fl\£). The 
continuous dependences of the solution none and T are given by theorem 5.1 of [ ], 
while the continuous dependences on uq and / are part of theorem 4.1 . 

This completes the summary of the linear theory and we now turn to the quasi- 
linear problem. 

3. The Quasi-Linear Problem. In [ , section 6] we have already shown that the 
quasi-linear problem 

(c[u}(x),Du) = f[u](x) ■ £ 2 in n\S , 

w|an = "0/ (3-1) 
(c[u](x),X7T(x)) > 6- |VT(x)| , 

where the dependencies of c and f on u are of a general functional type, admits a 
solution in BV(O). We achieved this result as a consequence of the Schauder fixed 
point theorem applied to the operator U [v] , which denotes the solution operator to 
the family of linear problems 

(c[v](x),Du) = f[v](x) -C 1 in Q\£ , veXcBV(n) 

u\dn = u , (3.2) 
(c[v](x),VT(x)) > B-\X7T(x)\ . 

In addition, we have given an example of non-uniqueness which tells us that we 
need stronger requirements on c and f in order to obtain uniqueness. 

The requirements used here which will suffice to ensure uniqueness and beyond 
that continuous dependence are summarized in the following. A special ingredient 
is causality w.r.t. the functional argument v, that is the dependence on v is of Volterra 
type w.r.t. the notion of time induced by our time function T. 

Definition 3.1. (Functional causality) Let T(x) be the time of the point x S O. The 
set Ciji x \ := {z 6 O : T(z) < T(x)}, denotes the past of the point x (w.r.t. the time 
function T). Let J-(Ci) and Q(Ci) be function spaces defined on Ci and let f be an operator 
of type 

f : J"(n) g(Ci) , with /[•](*) : ^(H) -> R d , xeCl. 

We say that f is functionally causal (w.r.t. time T) if 

f[v}(x)=f[vln T{x) ](x). 

This definition says that the value f[v] (x) depends only those values v(z) which 
V attains on the past fij-^) . Beyond that the functional causality implies the following 
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domain restriction feature. Let A be in the range of T. Then, for x £ Cl\ = {z £ Q : 
T(z) < A} and v £ ^(n), the inclusion ^t(x) ^ implies 

/M(*)=/[i7.]lo TW ](*)=/[i;.loJ(*). 

Hence, the domain restriction (onto Ci\ ) by f\ : .F(flx) — > Q(Q\), fx[v] = /Mln» is 
well-defined. 

Requirement 3.2. (Transport fields) Regarding the functional argument, transport 
fields are maps of type 

c : L l {D.) -> C 1 (n\S) 2 , anffe c[ . ] (x) : L 1 (n) — > R 2 , 

and n;e assume them to satisfy: 

1. c is functionally causal. 

2. ¥ or fixed v £ L 1 (D) the function c[v] : 0\E — >■ R 2 zs a transport field according to 
requirement 2.5. 

3. Uniformity of the unit speed and causality (iv.r.t. x) condition: 
a) \c[v](x)\ = 1 for all x £ fl\E and for all v £ L a (n). 

W T/iere zs « uniform lower bound f> > smc/z f/zai 

jS < (c[v](x),N(x)) < 1 Vi6n\E and Vu£L x (n). 

c) Bof/z conditions hold for the one-sided limits of c[v] on the relatively open C 1 arcs Ejt 
<>/X 

4. Bounds and continuity: 

a) The map D x c : L 1 (D) — > C(D\S) 2x2 - the derivative of c[v] w.r.t. the variable x - 
is L 1 bounded by 

\\ D xc[v]\\ L i {n) < Mi Vr?£L 1 (n)- 
W f/ie bound on \D x c[v] (x) \,from requirement 2.5 3a), 

\d x c[v] (x) | < m e , Vxen\y £ 

/zo/ds uniformly for all v £ L 1 (Q). 
cj f/ie map c zs Lipschitz in the following manner: 

||cM-c[u7]||oo < Li ■ ||p-tp|| L i(n) ■ 

Requirement 3.3. (Right-hand sides) Regarding the functional argument, right- 
hand sides are maps of type 

f : L 1 {Vl) -> C 1 ^) , wzf/z /[ . ] (x) : L 1 (Q) — > ]R , 

and we assume them to satisfy: 

1. f is functionally causal. 

2. Bounds and continuity: 

a) The map f is bounded by 

||/MIIco<m 2 Vi?£L 1 (n). 
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b) The map V x f : L 1 (D) — > C(fi) 2 - the derivative of f[v] w.r.t. the variable x - is 
bounded by 

||V*/[o]||oo <M 3 VoeiHn)- 

c) the map f is Lipschitz in the following manner: 

II/W-/HIU < U. - Ik-w|li,i(n) ■ 

Finally, we define, as in [ ], the subsets of BV(dCi) and BV(Ci) with which we 
will work later on. 

Definition 3.4. Let M\, M.%, M3 be the bounds from the requirements stated above. 

a) We denote by 

«8 = Q5(3n) := {v G BV(dn) : \\v || L oo (an) < M 4 , \Dv\ < M 5 } 

the set of boundary functions. 

b) Let the constants M* 6 1R and M** G IR be given by 

M * : =( M4 + /T^)- £2 ^' 

M 2 

+ J 

We set 



p—-. [ M 1 + \\DN\\ Vin]) 



X = X(n) := {v £ BV(Cl) : |M| L i (n) < , |Dr»|(n) < M„} . 

When speaking about transport fields c and right-hand sides / in the following, 
we mean transport fields and right-hand sides according to requirements 3.2 and 3.3. 

In [ ] we have shown that the solution operator IT : X — > X is a self-mapping. 
Here, U inherits additionally the functional causality of c and /. 

Lemma 3.5. The solution operator U:I->X of problem (3.2) is functionally causal. 

Proof. Let v G X be arbitrary but fixed. Then, according to equation (2.4), U[v] (x) 
is given by 

U[v](x) = uo(rj[v](T (x),x)) + / f Q [v] o rj[v](r,x) dr , 



where the backward characteristic rj [v] ( . , x) is the solution of 

y' = -c [v]{y) , y(0)=x, 
while Co and fo are given by 



(c[v](x),VT (x)) ' JVI n ' (c[v]{x),VT (x)) 
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Clearly, eg and /g are the same way functionally causal as c and /. 

Since Tg has the same level-sets as T while merely the names of level lines have 
changed, we refer to Tg in order to denote the past C1 To m := {z £ O : Tg(z) < 

T (X)}. 

For every t £ ]0, Tg(x)[ we know that t][v](t,x) £ Cl Tg r x y The latter is a conse- 
quence of the uniform causality w.r.t. x (requirement 3.2 part 3b)). Hence, n[v](. ,x) 
depends merely on the restriction of the transport field Cg [v] onto the past ^t (x) • By 
the functional causality, c [v] | n T w.r.t. y depends only on c ■ ln T . Consequently, 
for f £ ]0, Tg(x)[ the expression w.r.t. v also depends only on v ■ In y 

Thus, by the representation of above, it is obvious that U is functionally 

causal: 

U[v}(x) = U[vt nToix) ](x). 

□ 

Using the domain restriction feature, as discussed above, we define X A : 

Definition 3.6. Let A be in the range ofT. We denote by X A = X(D A ) := {i>|n A : 
the domain-restricted version of X from definition 3.4. 

In the following, we consider the domain restricted quasi-linear problem on C1 A : 

{c[u]{x),Du) = f[u]{x) ■ C 1 in Q A , 

"lan = "o, (3-4) 
(c[u](x),VT(x)) >p.\VT(x)\, 

The question for existence has already been answered in [ , section 6]: by lemma 3.5 
the operator U : X — »■ X is functionally causal. Hence, the domain restricted version 
U\ : X\ — > X\ is well-defined and U\ [v] solves the domain restricted linear problem 

(c[v](x),Du) = f[v](x) ■ C 2 in D A , v £ X A 

«|an = wo, ( 3 - 5 ) 
(c[p](*),Vr(je)> >/3-|VT(x)| . 

Since U A : X\ — > X\ is continuous (w.r.t. the BV weak* topology) and X A is non- 
empty, convex, and compact (see [ , section 6] for both statements), we get a solution 
of (3.4) as a fixed point of U A by the Schauder-Tychonoff theorem. And vice versa, 
every fixed point u of the original operator U:I->I after restriction u | q belongs 
to X\ and is a fixed point of U A : X A — >■ X A . 

In the following, we show that for every A in the range of T the fixed point of U A 
is unique. 

4. Uniqueness of the Fixed Point. The subject of this section is to show that the 
solution of problem (3.1) is unique or rather that U has a unique fixed point. For this 
purpose we will show that, for any choice of < A < 1, the domain restricted oper- 
ator U\ : X A — > X A is Lipschitz. Moreover, we will see that !i A is in fact contractive 
for a suitable choice of A. The contractiveness will then be the key to the uniqueness 
of the fixed point u — U[u]. 

In order to estimate the difference U [v i] — U[v2\ we prepare by setting up a PDE 
which is satisfied by the difference. For the purpose of abbreviation, we set 

c x := c[vi\ , c 2 := c[v 2 ] , and f x := f[v{\ , f 2 := f[v 2 ] ■ 
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Let u\ and u 2 respectively denote the solutions of the two linear problems 

(ci(x),D«) = f x {x) ■ C 2 in fi A , «|an = "o,i/ 

and (c 2 {x),Du) = f 2 (x) ■ C 2 in Q A , "Ian = "0,2 • 

As in the proof of lemma 3.5 we refer to the transformed time To instead of T and 
denote by Q A the lower level-set of Tg 

n A := {z G n : T (z) < A} , Ag]0,1[. 

For the first considerations we use different boundary data. When setting — 
u 0,2 = u o later on, we will obtain the relations 

u l = U[o 1 ], u 2 = U[v 2 }. (4.1) 

Let w denote the difference w := u\—u 2 . After having subtracted the problems from 
each other, the difference w must satisfy the linear problem 

(ci(x),Dw) = (ft(x) -fz(x)) ■ C 2 - (ci(x) -c 2 (x),Du 2 ) in n A , 
w\ dn = w , 

with boundary data Wq = «o,i — «o,2- 

By the same argumentation as used in [ , section 5] we see that w is the unique 
solution of this PDE. But, in order to solve for w, we cannot directly apply the method 
of characteristics as in [ , section 4], since the right-hand side is not an absolutely 
continuous measure. 

Instead we approximate the right-hand side by absolutely continuous measures. 
Since u 2 G BV(Q), there exists a sequence (u 2/ „) ne t^ of C°°(n) functions which con- 
verges strictly to u 2/tl , i.e., 

\\ u 2 — M 2,«|lL 1 (n) ~* an d ||jD«2|(n) — |Du2, n |(n)| — > . 

For this statement see e.g. [3, theorem 3.9]. Moreover, we have Du 2r „ = V«2, n (x) ■ C 2 
and |Dw 2 ,„|(n) = || Vu 2 ,„ || L i(n)- 

Using such a sequence we obtain an approximate problem 

{ci(x),Dw) = (Ji(x) -f 2 (x) - (ci(x) -c 2 (x),X7u 2 ,n(x)}) ■ C 2 in n A , 
Mm = , 

with a sequence of solutions w„ which we can construct using the method of char- 
acteristics. By scaling PDE (4.2) by the factor 1/ (c%, V7b), the family of forward 
characteristics £( .,s) is then given by the I VP 

j/ = ci,o(y)/ y(o) = ?(s). 

Here, we set 

r — Cl ( ■- h ar >A ■- Cl f0 — h 
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(Note: if we were to be consistent, we would set f 2i o := f 2 / ( c 2/ VXo), which differs 
from ■) According to equation (2.7), we obtain w n in characteristic variables by 

t 

w„ o£(f,s) = 7*w (s) + J (/i, - f 2 ~ (cyo ~ c\, Vw 2/ «)) °£(t,s) dT . (4.3) 

o 

The consideration of the sequence w n will not be of any use if zv„ does not tend to 
w in an appropriate fashion. We will show the desired convergence in lemma 4.2. 
But first, we rewrite w n o £ in order to get a more convenient representation of the 
difference (w n — w) o £(£, s). 

Because the PDE for U\ has the same transport field C\, we write U\ in character- 
istic variables as 

f 

mio£(£,s) =7*« 0/ i+ J fi i0 o£(r,s) dT . 



And because u 2/ „ is smooth, we have 

f 

«2,n°£(f,s) - w 2 ,„ o 7(s) = / (c 1/0 , Vm 2 ,„) o£(t,s) . 



Together with zt>o = mq,i — «o,2/ the last two observations imply 

w n o £(f,s) =u t o g(t,s) - u 2/n og(t,s) + u 2/ „oj{s) - J*u 0/2 (s) 
t 

'[4,Vu 2 , n ) - f$) o£(t,s) dT . 



Subtracting w = U\ — u 2 finally, we end up with 

(w n - w)o%(t,s) = (u 2 - u 2; „) o£(f,s) + (7*m 2 ,h(s) - 7*u 0/2 (s)) 



+ j ((4V%)-^)o?(v) 



dT. 



(4.4) 



As a second step of preparation, we show that requirement 3.2 part 4 b) implies 
uniform bounds on the determinant of D£. 

Lemma 4.1. Let cbea transport field. For fixed v e L 1 (D) let £ [v] denote the general 
solution of the IVP 

y' = coM(y) # y(0) = 7(s), urtft co[p](y) = / c[p ]^ vr (y)) ' (45) 
Trier/, 

flj £[z>] : ]0,1[ X ]a, b[ — > Q\(S U E) zs a diffeomorphism, where S := {£[z>](£,a) : f 6 
] 0, 1 [} . Here, f/ze domain ofs 6 ] a, & [ refers to f/ze domain of the boundary parametrization 
7, zfMe the time domain t 6 ]0, 1[ refers to i/ze range of Tq. 
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b) for < A < 1 there are bounds k\ and K\ such that 

< k A < detD£[v](t,s) < K A V (f,s) G ]0,A[ x ]a,b[ . 

The bounds A: A and K A depend only on A, but not on v. Moreover, A: A decreases, while K\ 
increases monotonically with A. 

Proof. Statement a) is taken from [7, corollary 3.3], so only statement b) needs a 
proof. On the restricted domain O a there is - by requirement 3.2 part 4 b) - a uniform 
bound on the derivative of c[v] 

\d x c[v](x)\ < M\ , Vxen A , VeeL^n). 

A similar bound - which we also call M A - will hold for the derivative D x Cq [v] (x) of 
the transformed transport field c [v] . 

Because g[v] solves the IVP (4.5) and |co[c](y)| < 1 / niQ f>, we have a first estimate 

detD£M < \U\v\\-\*£\v\\ < MV]1 



P ■ m 

The derivative d s £[v ] is the solution of the variational equation below 

deftv]' = DMA ° £M ■ y(o) = t'(s) ■ 

For every t e ]0, A [ the point £(f,s) belongs to O a . Thus, we estimate 

( 

\d s Z[v\\{t,s) < ||7'||co + J M A -\d s Z[v]\(T,s)dT, 


and an application of Gronwall's lemma (see e.g. [ ]) yields 

\d s Z[v]\(t,s) < ||7'||coexp(A-M A ) . 
Hence, we obtain the upper bound K,\ 

detD g [ P ]M< "^^J o A - M ^ =:X A 

on]0,A[ x )a,b[. 

For the lower bound we consider the inverse map for x £ O a \S. Ac- 

cording to [ , corollary 3.3], the inverse map is given by 

= (T {x),s[v](x)) , with s[v]{x) = ^ X {n[v} (T (x),x)) . 
Therein, rj [v] ( . ,x) denotes the backward characteristics given as the solution of 

y' = -c M(y) y(0) = xen A \s. 

Now, the determinant of D. Y £[c] _1 is bounded by 

dAD^[v]-\x) < |VT (x)| ■ \W x s[v](x)\ 



14 



T. MARZ 



and V x s [v] (x) T is given by the expression 

(j-'Yivf- (d,v[v}(T (x) r x) -VT (x) T + D xV [v](t,x)\ t=To{x) 

The derivative |DxJ/[f](i,3c)| can be estimated in the same way as |9 s £(f,s)|. And 
because x G fix, we obtain 

\D x t][v](t,x)\ < exp(A-M A ) . 
Using the latter result to establish a bound on | V^sfl?] (x)\, we end up with 



detDM-H*) < " VT °" L : ( ?;\| ■ ( 



exp(A-M A ) , 



for x G Oa\S. Finally we define l/fc A equal to the right-hand side of the last inequal- 
ity and get the desired lower bound 

detDg[ P ](t / 5) =det G Dg[p] ^^" 1 ) = detD *fW- 1 WI-ffM 

for (£,s) 6 ]0, A[ x ]a, since in this case £(f,s) 6 H A \S. 

Both bounds do not depend on the choice of V G L 1 (fi) and the monotonicity 
properties of fc A and K A as functions of A are obvious. □ 

In the next lemma we turn to the approximation of w by w n . 

Lemma 4.2. Let w and w n be as defined in the preparatory step above. Interpret the 
L 1 (fi x) -functions w and w„ as absolutely continuous measures w{x) ■ C 1 and w n {x) ■ £ 2 
on D A . Then, the sequence of measures w„(x) ■ C 2 weakly* converges to w(x) ■ C 1 : 

w n (x) ■ C 1 w[x) ■ C 2 , as n — > oo . 



Proof. Let cp G Co (O a ) be a test function. By changing variables it follows that 

b A 

J (w n — w)(x) ■ q>(x) dx = J J ((w n — w) ■ cp) o £(f, s) ■ detD£(f,s) dt ds . 

n A a o 

We use the representation of (w n — w) o £ according to equation (4.4) and study the 
convergence of the three summands in equation (4.4) separately. The first summand 
is estimated by 



b A 

J A(M2,»-w 2 )-<p)°£M-detD£(f,s) dt ds 

a 



— 1 1 ^2 n — ^2 1 



Li(n A ) 



and the right hand side tends to zero, because the sequence «2, n strictly tends to «2 
in BV(Cl). For the second summand we write 



J J (7*"2,«(s) - 7*"0,2(s)) ■ <po£(t,s) ■ detD£(£,s) dt ds 

a 

b 

/(7*"2,„(s)-7*"0,2(s)H7'(s)| 



7 A <p°g(f,s) -detDgfcs) dt 
IY0OI 



ds 
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Let k\ and K\ be the bounds on the determinant as in lemma 4.1. By the definition 
of k^, we have 

W(s)\ 1 1 

k K < detD£(0,s) < ip^i < 



6-m l?'( s )l fi-m-k\ ' 

And consequently, 

A 

/ fo^s) -detD£(£,s) df 
By the last result, we further estimate: 

b 

< A ■ ■ ||<p ||oo ■ / h*(u 2 ,n ~ "0,2)(s)| ■ |?'(s)| ds 







■ mo ■ 












■ m ■ 






Ka 






■ m ■ 





< A 

V>-m a -k x 



J I («2,« - u 2 )\do(x) | dH l (x) 



'lloo ■ || («2,w - u 2)|3nllLi(3n,-Hi) 
In the last factor we apply the trace operator for BV-f unctions 



I an 



which, according to [3, theorem 3.88], is continuous w.r.t. the strict topology on 
BV(Q). Hence, the factor ||(w2,« — ^lanll^fan/H 1 ) a ^ so ten ds to zero as n tends to 
infinity. 

Letip(f,s) := <po£(f,s) • detD£(f,s). Then, by changing the order of integration, 
we get for the third summand 



b A f 

' '„0 



8 

6 A 



cl,Vu 2 , n )-f2 )o?(T,s) lp(t,s)dt\ dTds. 



a 



By the definition of i/> and since £ is a diffeomorphism, there is a continuous function 
A £ C(n A ) such that 



/zo£(t,s) = (detD^s))- 1 - j Jip(t,s) 



dt 



With h, we write the last integral as 

b A 

' 'A, Vu 2 ,„) - ft] o £(t, s) ■ ft o £(t, s) ■ det Df (t, s) dr ds 



a 



/ ((^V« 2 ,„)-/ 2 °) (x)-ft(x) dx 
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Next, we use the fact that 1/2 solves the PDE (c2(x),Dti2) = fo(x) ■ £ 2 to formulate 
the last result as 

= J (h(x) ■ c\{x), Vw 2 ,n(^)) dx- J (h(x)>c%(x),dDu 2 (x)^ 
Ox n A 

= J (q>(x) ,V U2 >n (x)) dx — J (<p(x),dDu2{x)} . 
n A n A 

In the second equation we have set <p(x) := h(x) ■ c9(#) as a new test function which 
belongs to C(O a ) 2 . 

Owing again to the strict convergence of «2,n to 1*2, we argue by [ , proposition 
3.15] that the last integral expression tends to zero as n — »■ 00. Summarizing the three 
steps above we obtain 

J (w n — w)(x) ■ q>(x) dx — >■ V cp G Co(n A ) , 



which means w n 

Based on the properties of the sequence zv n we obtain an estimate of the differ- 
ence XV — U\ — «2 in the next lemma. Later on, this estimate will account for the 
operator U\ to be Lipschitz. 

Lemma 4.3. The difference w = u\ — 112 satisfies 



\w 



| L i(n A+/ ,) ^ + • C A+; , ■ ||u ,i - UQ i2 \\ L i^ nM i) 
+ C A+h -£ 2 (n) ■ ((A + /z) ■ HA -/ 2 || L - ( o A )+fe - ll/i -/2lk-(n A+ *)) 
+ C A +/, ■ M** ■ f(A + /i) ■ || ci - c 2 || L ~(n A ) + || c x -cz\\v(q x+ 



Here, the factor C\ := ^. mQ \ A is an increasing function of A. 

Proof. We use the approximation of w by w n again. Because of the weak* conver- 
gence according to lemma 4.2 and because of the lower semi-continuity of the total 
variation w.r.t. the weak* convergence (for the semi-continuity of norms, e.g., see 
[ ]), we have 



\w\ 



w-L' 



(Oi) < liminf w n ■ C 1 (Q>) = lirninf Hwnlln^n 



and thus we can estimate ||H>n||ii(n A ) instead. Using the representation of w n by 
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equation (4.3), we obtain 

b A+h 



J J \w n \ °£(t,s) • det D£(t,s) dtds 

a 

</ / l7*(wo,i-wo,2)(s)|-detD^,s) dtds 



\ W n\\L\n k+h ) ~ 

o. 

i) A+fr 



a 

6 A+/7 f 

+ // J \f h o-f%\°aT,s)dT-detDat,s)dtds 

a 
6 A+fr f 

+ // /|(( Cl <°" c 2)' Vu 2,")|°£( T ' s ) dT-detDZ{t,s)dtds. 

a 

By arguing the same way as in the proof of lemma 4.2 for the first summand we get 

b A+h 

J J |7* ("0,1 - "o,2)(s)|-detD£(f,s) it cfe < (A + /z) ■ C x+h ■ ||«o,l - %2llii(an,«i) ■ 

a 

For the last summand, let g(r,s,t) : = (|c 1/0 — c\\ • |V«2,nl) ° ?( T / S ) ' detD£(f,s). 
Then, we estimate 

b A+h t 

J J y^((c L0 -c!?),V U 2,,,)|o£(T,s)dT-detD£(f,s) dtds 

a 

b A+h t b A+h / A+h \ 

~ I I J g ( T,S,tS> drdtds = J J J g{t,s,t)dt\ drds 



a a \ T 

6 A /A+/i \ 6 A+h / A+h 



~ J J J S(^s,t)dt\ dTds + J J J g{r,s,t)dt \ drds. 
« \0 / a A V A 

For the inner integrals, we have 

A+h A+h 

J g(T,s,t)dt=(\c lfl -4\-\Vu 2 ,n\)oZ(T,s) J detD£(t,s)dt 



< (\c lfi -c° 2 \ ■ \Vu 2>n \)oZ(T /S )-{\ + h-T)-K A+h 

< (\cyo -4\- \Vu 2 , n \) o£(t,s) • detDf(T,s) ■ (A + h - r) ■ 

K A+h 



In the next step we take away the scaling factor, which is in the transport fields and 
the right-hand sides of the PDE, by 1/ (c x , VT ) < 1/ (m ■ |S): 



< (jd -c 2 | ■ |Vu 2 ,„|) o?(t,s) -detD£(T,s) ■ (A + Zt-T) 



&A+h 



p-m - k x +h 

(|c : - c 2 | ■ |Vu 2/W |) °£(t,s) -detD£(T,s) ■ (A + /i - t) ■ C A+h . 
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By the last result we infer on the one hand that 

b A A+h 

II I g(yT,S,t ^ dtdTds <( x + h )- c x+h J \ci~ c 2 \(x) ■\Vu 2 ,„\(x) dx 
a o o n A 

< (A + /i) -Cx+hWci -c 2 ||i»(n A ) ■ ll VM 2,n|lii(n A ) ' 

and on the other hand that 

b A+h A+h 

l l l g ( T,S '^ dt dT dS ~ k ' CX+H J l C l ~ C 2l( X ) ' l Vu 2,n|(^) ^ 

a A A Cl K+h 

< h ■ C A+h \\c 1 - C2|| L oo (nA+A) ■ ||VM2,„|| L i(n A+(l ) ■ 

Finally, for the summand 

b A+h t 

\fl,0 ~fi\° £(t,s) dT ■ det Df(f,s) dt ds , 



a 

we need to perform the same steps with 

g(r,s,t) := (|/i /0 -/ 2 °| ■ ln A+ J ^( T ' s ) ■ det D^(f,s) , 

and end up with the same estimates, but \\c\ — c 2||L°°(n A+;i ) has to be replaced by 
\\fl-fzh<"(n A+h ) an d l|Vw2,n|| L i(n A+ „) has to be replaced by £ 2 (n A+/l ). 

Summarizing the last considerations we have an estimate for \\w n || i^i (o A+/ ) by 

\\ l °n\\L\Cl A+h ) < ( A + ^)C A +/|- II "0,1 - "0,2|| L l(9aWl) 

+ C A+h £ 2 {FL A+h ) ((A + /i)- ||/i-/2|| L ~(n A )+ ft - \\h -/2|| L »(n A+/l )) 
+ C A+/! || Vm 2 ,«|Il1(o a+/! ) (( A + k) ■ \\ci-c 2 \\ L ^^)+h- \\cj 

~ c 2 |Il°°(O a+;! ) J ■ 

Because u 2/ „ strictly tends to u 2 , we have 

l|Vw 2 ,„|| L i(n A ) -> |DM 2 |(n A ) . 
Hence, going over to the lim inf and plugging in the estimates 

|DM 2 |(n A ) < M„ and C 2 {D. A+h ) < £ 2 (Q) 

finally yields the assertion. 

Because by lemma 4.1 we know that K A increases while k A decreases with A, it is 
clear that C A increases with A. □ 

As a corollary, we obtain that U A is Lipschitz. 

Corollary 4.4. For any choice of < A < 1 the operator U A : £ A — >■ X A is 
l}-Lvpschitz 

||irN-U[t»2]||z,i(n A ) < A ' K A- l^l-^|| L i(n A ) ■ 
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Here, k a is defined by 

K A :=C A -(L 2 -£ 2 (Ci) + L 1 -M^ , 
and is an increasing function of X. 

Proof. Let V\, v 2 G X A . When we consider the operator U or !i A , we always 
use the same boundary data uq £ *B. Hence, as mentioned in equation (4.1) of the 
preparatory part at the beginning of this section, we have w = u\ — u 2 = U\[vi] — 
IT A [v 2 ], since Mo,i — w 0,2 = u o- By using lemma 4.3 with h = we see that 

\\U[vx] - U[v 2 ]\\ L i {nx) < A-C A ■ (r 2 (n) ■ H/j -/ 2 || L ~(n A ) ■ Iki - c 2|| L ~(n A ))- 

For the differences f\ — fi and C\ — c 2 we use now the functional causality and the 
Lipschitz conditions, which we require. That is 

ll/i-/2||jr°°(n A ) = ll/N -/Mli°°(n A ) = II (/N -/N)|n A l|oo 

= ll(/h-lnJ-/b2-lnJ)|n A ||oo< HM-lnJ-Z^-lnJIU 
< L 2 ■ ||t?i ■ ln A - v 2 ■ ln A |lii(n) = L 2 • ||oi - »2lLi(n A ) 

and analogously 

Iki - c2|| L ~(n A ) < L i ■ \\ v i ~ p 2lli,i(n A ) ■ 
Putting everything together, it follows that 

|| U^] - !iMI Ll( n A) < A ■ C A ■ (L 2 ■ £ 2 (Q) + Lj ■ M„) ■ ||n - ^|| L i(n A ) 

< A-k a ■ \\vi -»2|| L i(n A ) • 

Finally, jc a increases with A, since C A does so, too. □ 

Now that we have brought together all ingredients we are able to show the 
uniqueness of the fixed point. 

Theorem 4.5. (Uniqueness) The solution operator U : X — >■ X of the (non-restricted) 
original problem (3.2) has a unique fixed point u £ X, u = U[u]. 

Proof. First, we show that, for any choice of < A < 1, the domain-restricted 
operator U\ : X\ — > X\ has a unique fixed point. In order to do so we decompose 
Q A into finitely many stripes Q/j^jy), of "thickness" h > 0. For such stripes, we use 
the notation 

n A+M = n A+/; \n A = {z e n : a < r (z) < A + fe} . 

Let the step size h be such that h < 1/k a , and let L = |_A / /xj 6 N be the number of 
steps. Then, 

L-l 

= U ft(l+l)h,lh U ^\,Lh • 
/=0 

For the first step, consider the operator : ■£/, — > ■£/, on CI;, = CI/, - By corollary 4.4 
and the choice of h we have a contraction 

\\U h [vi] - U h [v2\\\ L i {Cth) <h-x h - \\v 1 -v 2 \\ L i {nh) <h-K X - \\vx-v 2 \\ h i { ci h) . 
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If now u\ — U[u\] and u 2 = U[u2\ are two fixed points, we have, after domain- 
restriction onto Cl h , 

II "l ~ u i\\v-(n h ) < h ' K \ - - "2|Li(n ft ) 

and consequently 

0< (l-fe-K A )||ui -U2\\ L i { a h) <0. 

Hence, all fixed points coincide on the stripe O/, . 

Next, we perform an inductive step. Assume that all fixed points coincide on 
Q;/,, we show that they must also coincide on On+ijh = n;/ i+ ;,. Let u\ = U[u{\ and 
u% — U[u2\ be two fixed points again. With w = u\ — u 2 , by lemma 4.3 we know that 

ll"i-"2|| L i(n (m)(! ) <C(z+i)fc- 

(c 2 (Cl) ■ ((/ + l)h ■ \\f[ Ul ] - f[u 2 ] || t . (0a) + /z ■ - /[« 2 ] k-(n (m)fc) ) 

+ M„ ■ ((Z + l)fc. ||c[ui] -c[M 2 ]|k»(n tt ) +h- \\c[ui] -c[u 2 ] || L»(n (m)fc )) 
Because «i and w 2 coincide on fi^, we have 

ll/N -/MllL~(n tt ) =0 and ll c l"i] - c ["2]||L~(n„,) =0, 
and thus, the estimate reduces to 

ll"i -"2|| L i(n (m) „) < C(i+i)h-h- (£ 2 (n) ■ \\f[ui] - f[u 2 ]\\ L ^n {m)h ) 

+ M„ ■ ||c[tt X ] -c[u 2 ] || L co (n(I+1)(i) ) . 

By using the Lipschitz conditions on c and / again we have 

ll"i " "2|| L i(n (m)4 ) < ^(i+i)k ■ ll«i - "2|| L i(n ((+1)fc ) 

< toc A - - M2|ll.l(n a+1)Ji ) • 

More precisely, because we have assumed that U\ and w 2 coincide on D H „ the latter 
inequality in fact means 

ll"i " "2ll Ll( fi (w ) < hK A • ||m - "2|| Ll( n (m)WA) • 

By the contractiveness, < 1, we see that ||mi — W2||L 1 (n (/+1)/ u ) = and so me 
fixed points also coincide on the next stripe ,//,■ 

For the last stripe we have to adapt the step size to h = A — Lh. But, since h < h, 
the same argumentation applies. 

As claimed before, the domain restricted operator U\ : X A — >■ X 7 \ has a unique 
fixed point for any choice of < A < 1. Now the last step: assume by contradiction 
that the non-restricted operator U : X — >■ X has two different fixed points, U\ and m 2 . 
Therefor, u\ and m 2 must differ on a subset W c CI with £ 2 (W) 7^ 0. Because the 
stop set £ has Lebesgue measure zero, £ 2 (£) = 0, we can choose < A < 1 so close 
to 1 that £ 2 (n A nff)/0. Thus, we have 

||"i - "2||i,i(n A ) = II "1 ~ M 2|lLi(n A nw) 7^ • 
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But, because «i |n A and «2|n A are fixed points of the domain restricted operator iI A : 
X A — > £ A , we also have 

IN - W2|| L i(n A ) = 

by the previous uniqueness proof. A contradiction. □ 

5. Continuous Dependence of the Fixed Point. In this section, we will show 
that the unique fixed point depends L 1 -continuously on the following data: the 
transport field, the right-hand side, and the boundary data. We consider two lin- 
ear problems: 

(c[v](x),Du) =f[v]{x) -C 1 in Q\E, 
"Ian = "o / 

and 

(c[v]{x),Du) = f[v](x) -C 1 in Q\E, 
"Ian = "o / 

where we assume that for both problems the same domain Cl, the same stop set E, 
and the same time function T (with transformed version To) are specified. 

Moreover, we assume that c and c both satisfy the requirement 3.2 with the same 
bounds, and that f and f both satisfy the requirement 3.3 with the same bounds. 
Finally, we assume uq G 03 and Uq G 03. By the latter assumption we are sure that we 
obtain two solution operators 

U:X^X, v^U[v], 
U:X^X, v^U[v], 

which respectively correspond to the two linear problems above and possess the 
same domain and range X, which depends on all those bounds. 

We view c and / as perturbed versions of c and /. In order to measure the 
perturbation we introduce the following norm: 

Definition 5.1. For maps of type g : L 1 (n) — s- Cb(Cl\L) d or of type g : L^O) — s- 
C{Cl) d , d G N, we define the norm 

\\g\\o-= sup HsMlloo. 

Theorem 5.2. (Continuous dependence) Consider the two solution operators U and U 
as described above. Let u = U[u] and u = U[il] be the unique fixed points of these operators. 

Then, for every e > 0, one canfind 3 > such that \\u — u\\ l i^q\ < e, whenever 
(ll"o- "o|| L i(3fW) + £ 2 (°) • 11/ -/Ho + M**-||c-c|| ) < S. 
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Proof. Let v,v £ X be arbitrary but fixed, and set «i := U[c] and u 2 := As 
before, we set up a PDE for the difference w := U\ — u 2 on the restricted domain Ci A , 
< A < 1 by 

(c[v]{x),Dw) = (f[v]-f[v])(x)-C 2 -{(c[v}-c[v})(x),Du 2 ) in n A , 
^lan = "o - "o • 

Again, we choose a sequence u 2/ „ 6 C°°(n) which strictly converges to u 2 in BV(C1). 
And again, the sequence w n of solutions to the approximate PDE - which has the 
absolutely continuous measure Vu 2/H (x) ■ £ 2 instead of D«2 - converges weakly* to 
w. 

In order to proceed as in lemma 4.3 we supplement the right-hand side of the 
approximate PDE to 

(c[v](x),Dw n ) = {(f[v] -f[v}) - (c[v]-c[v},Vu 2 , n )) (x) ■ C 2 

+ am -/>]) - m - m^^)) w ■ & ■ 

For the first summand of the new right-hand side we will apply the steps from the 
proof of lemma 4.3. For the second summand we operate in a simpler way. Let 
£ = f [v] be the characteristics corresponding to the field Cq[v] and let 

_ \(m-m)-m-m,vu 2 , n )\ 

g (cM,VT ) 
As in the proof of lemma 4.3 we have to estimate 

b A+h t 

go £(t,s) drdetDf(f,s) dt ds < ... . 

a 

After having changed the order of integration and having estimated the determinant, 
we obtain 

b A+h 

... < ( A + fc).£A+ft. / / go£(T,s)detDf(T, S )dTds 

K A+h J J 
a 

<(A + h)-C A+h - J \(f[v]-f[v])-(c[v]-c[v},Vu 2 , n )\(x)dx 

<(A + h)-C A+h (£ 2 (fi) ■ ||/[0] "/>]||oo + ||Vu2,„|| L i (n) ■ \\c[v] -WW*) 
<(A + h)-C A+h (£ 2 (n) ■ ||/-/||o + ||Vu 2 ,„|| Ll(n) ■ ||c-c||o) ■ 
Putting both estimates together and then going over to the lim inf, we end up with 

IM| L i (nA+(;) < (A + h) ■ c A+h ■ IK - «olln(aa«i) 
+ C A+; , ■ £ 2 (n) ((A + h) ■ \\f[v] - f[v] || L » (nA) + h ■ \\f[v] - f[v] || L ~ ( n A+(I )) 
+ C A+/l -M„ ((A + /z) ■ ||c[z7] — c[0] || L =«(o A ) +h- IkM -c[z)]|| L co (nA+(i) ) 
+ (A + h) ■ C A+h (£ 2 (Q) ■ ||/ - /||o + M„ ■ ||c - c|| ) . 
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Now we can show the continuous dependence in the domain restricted situation. 
Fix < A < 1, choose a step size < h < 1/k a and let L = [X/h\ 6 N be the 
number of steps. Furthermore, let 



u o - uoWiMdriM 1 ) + £ ( n ) ' 11/ -/Ho + M ** ■ lk-c||oJ < S, 
for some S > 0. Let / 6 No, I < L. With the result above we estimate on the set 

IIHlLi(n (m) „) < A-C A -<5 

+ C A ■ £ 2 (n) ■ (A ■ \\f[v] - f[v] || L . (ntt) + ft • \\f[v] - f[v) \\ L -(n mv ,)) 

+ C A -M„ ■ (a - \\c[v]-c[v}\\ L co {Ctih) +h- \\c[v] -c[v]\\ L c° {n(l+i)h] 

By using the functional causality, the Lipschitz conditions on c and /, and the defini- 
tion of jc a from corollary 4.4 we obtain 

IM| L i (fW) < A-C a -J + Ak a - \\v-v\\ L i {Cllh) +hK k - \\v-v\\ L i {C i {l+l)h )- 

Let 5 = A • C A • 5. Now, we plug in the two fixed points u and u, i.e., we set 
u\—v = u and ui = v = u, 

II" " "IlH^W) < ^ + A ■ k a ■ || U - u|| Ll(nft) +fc • k a • ||u - w|| L x ( n (m)fe ) ■ 

We define the error e; on the set Q/j, by e; := ||u — m||^i(q (( y Then, by our choice of 
h, the last estimate yields the error recursion 

< (1 - h%x) ■ e/+i < S + A ■ k a ■ e, , 

which leads to 



i 

e/+i < £ a* • i with 01 



5 , . A ■ k a 
K — with a := . 

HKi 1 — rtKl 



In summary, we get 

/l-a L+1 A-C A \ r 

and the continuous dependence for the domain restricted case is obvious. 
Final step: let £ > 0. For the full domain fi we choose A so close to 1 that 

e 

II" " u \\V-{Ci\Ci k ) ^ 2 ' 
In dependence of this A we find h and L. What remains to do is to require 

l-a L+1 A-C» \ _1 £ 



1 - a. 1 — hx\ J 2' 
then, we get ||w — «||jifn) < £, whenever 

(ll"o — "o|| L i(afW) + £2 ( n ) ■ 11/ -/Ho + Mw||c-c|| ) 



< J. 
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6. Discussion. We have shown in this paper that the quasi-linear problem (3.1) 
has a unique solution which depends continuously on the coefficients of the PDE and 
the boundary data. The special ingredients here have been the functional causality 
and the Lipschitz continuity of the coefficients. 

The good behavior of the time function T and the transport field c close to and 
on the stop set E have not been used in sections 4 and 5. So, in the theory presented 
we could consider solutions which belong only to B V(Q\E). 

In [7], the mentioned good behavior of T and c close to and on E accounted 
for closing the gap E: for the linear case we obtained that the solution belongs to 
BV (fi) and not only B V(n\E), see [ , theorem 4.1]. This feature was crucial to proof 
the existence of a solution to the quasi-linear case (even for non-causal functional 
dependence) by the Schauder fixed point theorem because X C BV(fi) is now com- 
pact. Moreover, the fixed point solutions belong themselves to BV(d). So here, we 
have again the closing the gap feature and additionally uniqueness together with 
continuous dependence. 

The advantage of closing the gap here is, that we can allow for more general time 
functions. For example we can use one with levels as illustrated in figure 6.1. Here 




Fig. 6.1. Three levels of a time function T with a locally maximal level in the middle (dashed green line). The 
white area is the domain O with its boundary, which is the start level T = 0, in red. The solid green and the dashed 
green lines together are a saddle level (the terminal nodes of the dashed green line behave like saddle nodes). The 
dashed green line alone is an intermediate stop set while the blue lines are the maximal level of T , i.e., the final stop 
set which is disjoint. 

in a first step one solves from the start level to the saddle level. With our closing the 
gap feature we can even reach the saddle level and close the first gap which is given 
by the intermediate stop set. Finally, one solves the remaining two problems and 
gets a unique global solution in fi. Without closing the gap, one could only produce 
a local solution which is defined on the region between the start and the saddle level. 
For more on this topic see [ , chapter 5]. 
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